Helium-4 Luttinger liquids in nanopores 
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We study the low temperature properties of a 4 He fluid confined in nanopores, using large-scale 
Quantum Monte Carlo simulations with realistic He-He and He-pore interactions. In the narrow- 
pore limit, the system can be described by the quantum hydrodynamic theory known as Luttinger 
liquid theory with a large Luttinger parameter, corresponding to the dominance of solid tendencies 
and strong susceptibility to pinning by a periodic or random potential from the pore walls. On 
the other hand, for wider pores, the central region appears to behave like a Luttinger liquid with a 
smaller Luttinger parameter, and may be protected from pinning by the wall potential, offering the 
possibility of experimental detection of a Luttinger liquid. 



While the ability to flow through narrow pores is what 
gives a superfluid its name, a strictly one-dimensional 
(ID) Galilean invariant system of bosons with short range 
interactions cannot exist in any true ordered (superfluid 
or solid) phase. Instead, it will be in a type of quasi- 
ordered phase known as a Luttinger liquid (LL), featuring 
correlations that decay as powers of spatial separation, 
even at zero temperature (T) [lj. Ongoing experiments 
on 4 He in nanopores, have recently moved towards pore 
radii in the nanometer regime, thus offering the exciting 
possibility of probing the LL phase. Historically, confinc- 
ment was achieved through quasi-lD cavernous networks 
in porous glasses [2j and more recently from folded sheets 
of mesoporous materials Q . A substantially different ap- 
proach has been taken by Gervais and collaborators [J], 
who are studying the flow of 4 He inside nanopores of cus- 
tom radii by carving a hole through a S13N4 membranc 
using a transmission electron beam. Although the exper- 
iments of Ref. [IJ have thus far focused on flow proper- 
ties of helium in the gas phase, it is intriguing to ponder 
the equilibrium system of helium atoms inside the pore 
as the temperature is reduced below the bulk superfluid 
transition temperature T\ ~ 2.17 K. If the pore radius 
is sufficiently small, it ought to be possible to observe a 
crossover to strictly ID behavior. 

Significant progress has been made on the theoreti- 
cal understanding of helium confined inside carbon nan- 
otubes Ø-Q or smooth nanopores [1,13 . A complex phase 
diagram has been predicted, containing states where he- 
lium atoms occupy only the central region of the cylinder 
(for narrow tubes) and those consisting of one or more 
cylindrical shells (for wider tubes). A realistic micro- 
scopic description of an assembly of 4 He atoms confined 
inside a single nanopore can be achieved by making use 
of the accepted Aziz [ll| pair potential to describe the 
interaction between two helium atoms. The pore can be 
modeled as a long cylindrical cavity carved inside a con- 
tinuous medium. The potential energy of interaction of a 
single helium atom in the pore can be obtained by inte- 
grating a Lennard- Jones pair potential, in the same way 



as it is done for smooth pianar substrates [12[. 

In an attempt to investigate theoretically the physical 
issues addressed by the planned experiments, in this work 
we made use of such a model to reproduce as closely as 
possible the proposed experimental geometry. An even 
more realistic simulation may require including the ef- 
fects of a periodic substrate, disorder in the pore walls 
and flow conditions, but understanding the equilibrium 
state with clean smooth walls is an important first step. 

We computed low-temperature thermodynamics of 
fluid 4 Hc inside nanopores by means of Quantum Monte 
Carlo (QMC) simulations, based on the continuous-space 
Worm Algorithm (WA) . This methodology affords an es- 
sentially exact estimation of many physical observables, 
for systems of thousands of quantum particles at low tem- 
perature, in the grand canonical ensemble H, 14 1. The 
possibility of simulating large systems is crucial, as LL 
behavior can only be detected on sufficiently long length 
scales. We chose Lennard- Jones parameters for the atom- 
wall potential appropriate for SiaN4 [Hj], as well as glass 
[l6| , and found no qualitative change in the basic physical 
results. All results shown here are for the S13N4 system. 

We considered nanopores of length L up to 128 Å and 
radii R between 2.5 and 12.0 Å, using periodic boundary 
conditions along the axis. The illustrative results shown 
here are all for £=100 Å; we observed that physical re- 
sults are insensitive to the length of the nanopore, pro- 
vided that L > 8R. Simulations were carried out in the 
grand canonical ensemble, for temperatures in the range 
T = 0.5 — 2.0 K, at a chemical potential fi = —7.2 K, 
which corresponds to saturated vapor pressure in the 3D 
reservoir used in the experiments. In addition to restrict- 
ing T < T\ ~ 2.17 K, all temperatures considered are 
small with respect to the average kinetic energy per par- 
ticle, ensuring we are in a low-energy regime. A typi- 
cal simulation employs an average number of 4 He atoms 
around 1000. 

Fig. Q] shows the computed density p(r) of 4 He atoms 
as a function of radial distance r from the tube center, 
at T = 0.5 K. Results are shown for pores of diffcr- 
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FIG. 1. (Color online) QMC data for the radial density 
of helium atoms for cylindrical nanopores with radii R = 
2.9, 4.0, 10.0, 12.0 Å showing the characteristic shell structure 
due to the interaction of bulk atoms with the walls of the 
channel. Error bars are smaller than the line thickness. 



ent radii. The density is uniform along the axis of the 
nanopore (z), and also as a function of angle. As ex- 
pected, p(r) — > in the vicinity of the wall, and features 
a large peak approximately 2.5 Å from its surface, due 
to the short-distance repulsion and pronounced minimum 
near the surface. One or more peaks in the radial density 
are observed, corresponding to a possible inner cylindri- 
cal region of high linear density surrounded by cylindrical 
shells of atoms separated by a distance set by the attrac- 
tive well of the Aziz potential. This phenomenon is very 
robust and always occurs for sufficiently large R. As R 
increases, more peaks occur, with those near the center 
of the pore gradually evolving to a constant density in- 
dicating that dimensional crossover to the 3D limit has 
occurred. The presence of an inner cylinder (IC) near 
the axis of the tube depends on the precise choice of R, 
and for the particular form of the helium-pore potential 
used here, occurs for R ~ 3n where n > 1 is an inte- 
ger. Cylindrical shells in a nanopore are the analogue 

illrzi- 



of pianar layers in a He-4 thin film [lg, |T7| . They arise 
from the suppression of quantum fiuctuations near the 
container wall, where atoms are localized fl6j . As men- 
tioned above, such layering has been predicted previously 
in nanopores using variational and approximate Density 
Functional Theory [1, 0|, as well as QMC at very low 
helium density |10|. Here, we observe the formation of 
up to four such concentric shells at p = —7.2 K, owing 
to the relatively large size of the systems simulatcd. 

An attempt to develop a theoretical understanding of 
the low energy properties of the helium atoms in the pore 
can be made through Luttinger liquid theory which, 
in strictly ID, provides a universal description of inter- 
acting fermions or bosons via linear quantum hydrody- 
namics. This is accomplished in terms of two bosonic 
helds, 8(x) and 4>(x), representing the density and phase 
oscillations of a second quantized particle held operator 
i)^{x) ~ po + d x 9(x)e l<l> ^ related to the /x-dependent 
density of particles at T = in the thermodynamic limit 



via po = (p(x)) = {ip^ {x)4>{x)) . The quadratic Hamilto- 
nian describing these helds is given by 
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where we have set H = &b = 1 and v and K are related to 
the microscopic details of the underlying first principles 
many-body Hamiltonian upon which the simulations are 
based. The velocity v describes the linear dispersion of 
low energy phonon modes, while the Luttinger parameter 
K can be tuned to initiate a T = crossover between a 
state with solid order at K = oo to onc with inhnitc range 
supcrfluid correlations at K = 0. The quadratic nature of 
Eq. ([1]) allows for the calculation of all correlation func- 
tions and thermodynamic properties in terms of v and 
K. The continuous-space WA has recently been success- 
fully employed 18[ to test the effectiveness of Eq. ([1]) in 
describing the low energy quantum dynamics of a system 
of interacting bosons in the ID continuum with contact 
interactions. The methodology of Ref. [l8[ is used here 
to study the density-density or pair correlation function 
(PCF), which can be derived for a LL to be: 
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with r/(iz) the Dcdckind eta function and 9i(y, z) the first 
Elliptical Thcta function. A is a non-universal constant 
dependent on the short-distance properties of the system. 
As LT/v oo, Eq. $2$ simplifics to (p(x)p(O)) -> pi - 
l/(2Tr 2 Kx 2 ) + Acos(2Trp x)/x 2 / K Ø. 

For R = 2.9 Å, Fig. [TJ shows that helium atoms are 
confined to the center of the pore (inner cylinder) with a 
radial density that is effectively zero by r ~ 1.25 Å; we 
thus expect ID behavior in this case. The main panel 
of Fig. [2] shows the computed axial PCF for tempera- 
tures ranging from 0.5 - 1.25 K. The results suggest that 
helium inside the narrow (R = 2.9 Å) nanopore is m a 
quasi-solid phase, with slowly decaying correlations only 
minimally dependent on temperature. At low temper- 
ature, the average effective ID density (po = N/L) is 
close to r^ 1 where ta — 2.97 Å is the minimum of the 
helium-helium interaction potential. In this ID limit, we 
expect the PCF to be well described by LL theory and 
can perform a fit of the QMC data to Eq. (J2|). The He- 
pore interaction is independent of the axial coordinatc, 
and Galilean invariance restricts the ratio v/K = irpo/m 
The finite size and temperature scaling behavior of 
all thermodynamic quantities is predicted by LL theory. 
Thus, a single fit of any quantity, performed for particular 
values of T and L, can be used to generate LL predic- 
tions throughout the region of T — L parameter space 
where the universal hydrodynamic theory of Eq. ([T]) ap- 
plies. The results of this fitting procedure are shown as 
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FIG. 2. (Color online) QMC data (symbols) and a fit to 
Eq. @ (lines) for the pair correlation function along the axis 
of a nanopore with L = 100 Å. The main panel shows strong 
oscillations for R = 2.9 Å, while the inset details rapid decay 
for R = 4.0 Å (here the line is a guide to the eye). Error bars 
are smaller than the symbol size and data in the main panel 
has been given a vertical T-dependent shift. 



FIG. 3. (Color online) QMC data (symbols) and a fit to 
Eq. @ (lines) for the inner cylinder pair correlation function 
along the axis of a nanopore with L — 100 Å and R = 12.0 Å 
for helium atoms with r < 1.75 Å. The inset shows the super- 
fluid fraction for the full nanopore as a function of tempera- 
ture. Error bars are smaller than the symbol size and data in 
the main panel has been given a vertical T-dependent shift. 



solid lines in the main panel of Fig. [2] we have deter- 
mined v 2 .g = 70(3) ÅK and K 2 . 9 = 6.0(2) where the 
subscript refers to the radius of the pore and the number 
in the bracket is the uncertainty in the final digit. This 
large Luttinger parameter corresponds to a strong ten- 
dency toward solid formation. The Fig. [2] inset shows the 
rapidly decaying PCF at T = 0.5 K for R = 4.0 Å. Here, 
no IC exists, and the helium atoms form a cylindrical 
shell (see Fig. Q]) with interaction energy minimizing he- 
lical density correlations at short distances that serve to 
wash out possible LL oscillations at longer length scales. 
We thus postpone a LL analysis of the R = 4 Å pore to 
future studies. In both these narrow pores, atoms expe- 
rience a large degree of localization, either along the axis 
or near the surface, and quantum-mechanical exchanges 
are strongly suppressed. 

The situation is markedly different for the R = 12.0 Å 
nanopore, as p(r) plotted in Fig.Q]indicates the presence 
of an axial IC of helium, surrounded by three cylindrical 
shells. Although the density between the shells is never 
strictly zero, in order to make direct comparisons with 
the R = 2.9 Å nanopore, we can measure the proper- 
ties of only those helium atoms inside the IC defined to 
include all atoms with a radius smaller than the loca- 
tion of the frrst minimum in p(r). We find that these IC 
atoms make the dominant non-background (oscillatory) 
contribution to the PCF and the assignment of a given 
indistinguishable helium atom to the IC is performed dy- 
namically whenever measurements are made. QMC data 
for the IC-PCF are shown as symbols in Fig.[3]for temper- 
atures ranging from 0.5 - 2.0 K. We observe persistent 



oscillations that decay much more rapidly than in the 
narrower pore, with an envelope that contracts as the 
temperature is reduced. This markedly liquid-like be- 
havior is consistent with the onset of a finite superfluid 
response inside the cylinder. The superfluid fraction can 
be determined through the winding of imaginary-time 
particle world lines, as they wrap around the system in 
the axial direction (with periodic boundary conditions) 
,19|. Results for the superfluid density computed in this 
way for the full R = 12 Å pore (L = 100 Å) are shown 
as an inset in Fig. [3] 20]. For comparison, the R — 2.9 Å 
nanopore of the same length shows no evidence of super- 
fluidity over the same temperature range. 

We may regard the R = 12 Å nanopore as a coupled 
multi-component LL, with cylindrical shells replacing the 
legs of previously studied ladders [21( . Similar to the con- 
clusion for a 2-leg bosonic ladder, we might expect that 
only one "center of mass mode" survives as a gapless 
degree of freedom in the low energy effective field the- 
ory, due to tunneling between the shells, manifest here as 
multi-particle quantum exchange cycles connecting them. 
Under this assumption, we have performed a fit of the 
lowest temperature QMC data in Fig.[3]to Eq. @. When 
considering the IC, we no longer have Galilean invariance, 
and must independently extract v and K from the data. 
We find V12 = 42(2) ÅK and K 12 = 1.3(1), which have 
been used to plot LL predictions (solid lines) that agree 
remarkably well with the QMC data up to T = 2.0 K. We 
stress that after the LL parameter and velocity have been 
determined for a single temperature, only one fit param- 
eter A remains to be determined at all other values of T. 
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The surprising robustness of the LL description of the IC 
can be tested further by using the values of V12 and Kyi 
determined from the PCF, to generate predictions for 
other quantities. We have compared the single-particle 
density matrix (tpi (x)ip(Q)) computed in the QMC with 
the expected result from LL theory (with no new fitting 
parameters) and find acceptable agreement at low T. 

It is natural to ask how the form of the confining poten- 
tial might affect the low energy behavior of the nanopore 
system. We have performed simulations of helium in the 
exact ID limit (R — 0) with /i = /j,^ — Vp(0, 2.9) where 
Vp(r, R) is the He-pore interaction for a pore of radius R. 
An analysis of the PCF yields values of v — 74(2) ÅK 
and Kq = 6.3(2) which are in relative agreement with 
those found for the R — 2.9 Å nanopore. As the radius is 
increased, the main consequences of the cylindrical shells 
of helium that surround the IC are to screen the inter- 
action with the wall and drastically alter the shape of 
the He-He potential. The latter effect can be quanti- 
fied by computing an effective ID interaction potential 
V m (z,R) = ^-Jd 2 Cjd 2 CV A (\f-^\)p(C) P (C) where 
Va is the Aziz pair potential, piD = J^ 2 Cp(C) t ne e f~ 
fective ID density, p(£) the radial number density and 
f = (£, 4>, z) in cylindrical polar coordinates. We find 
that Vid(z, 2.9) « Vx{z), consistent with our expecta- 
tion that this nanopore radius approximates the ID limit. 
However Vib(z, 12.0) exhibits a broad shallow minimum, 
shifted from ta with a depth that is 1/2 of that of the 
unscreened bulk potential Va. The softening of the effec- 
tive interaction potential decreases the energetic benefits 
of forming a solid with lattice spacing near ta and as a 
result the helium atoms delocalize into a smaller- K Lut- 
tinger liquid with a finite superfluid density. 

The value of A^.g found here in the ID limit is thirty 
times larger that that found for a screw dislocation of 
diameter 6.0 Å in solid 4 He [3]. While the difference in 
K- values may be due to the different chemical potential 
(p = 0.02 corresponding to the bulk melting point) or 
effects of the screw dislocation potential (which may vary 
greatly from Vp), it deserves further study. 

The numerical value of the Luttinger parameter K for 
pores of varying radius can provide important informa- 
tion on the sensitivity of the LL to perturbations com- 
ing from disorder or commensuration. In the ID case, a 
weak periodic substrate, commensurate with the density 
is only irrelevant for K < 1 /2 [lj while weak disorder is 
only irrelevant for K < 2/3 [22j |. We have found a value 
of K « 6 at saturated vapor pressure for the narrowest 
pores, indicating a strong tendency to form a solid, re- 
sulting from the shape of Va(?~)- This large value could 
in principle be experimentally tested along the lines of 
Ref. [1] as the formation of a quasi-solid would impede 
the flow of helium through the pore. Additionally, the 
results presented here may be relevant to the interpreta- 
tion of neutron scattering measurements of the momen- 



tum distribution of helium in porous media [23j, where 
we would expect the structure factor to exhibit a broad 
feature with intensity characterized by a power law de- 
pending on K at T = 0. 

As the radius of the pore is increased, the formation of 
shells near the pore wall may serve to screen the central 
region from the disorder and periodic modulation of the 
wall potential. Whether this effect, plus the related re- 
duction of the Luttinger parameter in the larger radius 
pores leads to a localization length longer than the tube 
length, is a crucial question for the experimental observ- 
ability of LL behavior. 
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